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Clarke and Pringle (2004) derived a proper viscosity formula in a rotating gas by applying 
mean free path theory. We study their argument in detail and show that their result can be 
derived with a much simpler calculational procedure and a physically clearer picture. 

§1. Introduction 

The subject of angular momentum transport is important in theories concerning 
accretion disks. The a model of Shakura and Sunyaev^^ assumes outward transporta- 
tion of angular momentum due to some viscous force. This viscosity is believed to 
originate from turbulence, because the molecular viscosity is far too small in this 
case. 

Because the viscosity formula for turbulence is not well known, generally ac- 
cepted practices have relied on the formula for molecular viscosity, with the viscosity 
coefficient assumed to be far larger than that for actual molecular viscosity. In the 
present paper we investigate the molecular viscosity formula, which can be unan- 
biguously defined, applicable to a rotating gas. The viscosity formula for shear flow 
is well known. The shear stress is proportional to the rate of strain. For a shear flow 
represented hy u = (0, U{x)), the x-y component of the shear stress is given by 

axy = -rj{dU/dx) , (1) 

where rj is the viscosity coefficient. 

We next consider a rotating flow represented hy u = (0, Rf2{R)) in cylindrical 
coordinates, where R and J? are the radial distance and angular velocity, respectively. 
In this case the R-(j) component of the shear stress is 



CTR^ = -rjR{d[2/dR). (2) 

Microscopically, the formula can be derived with mathematical rigor from the 
Boltzmann equation using the Chapman-Enskog expansion.^) It can also be derived 
more readily by applying the mean free path theory of the kinetic theory of gases 
heuristically.^^'^^ A simple application of the mean free path theory to a rotating gas, 
however, leads to the incorrect result oc —d{R?f2)/dR}^ The correct derivation 
starting from the Boltzmann equation naturally leads to the formula (Pj)-^-' '^^ It is 
a puzzle why the application of the mean free path theory, which seems physically 
plausible, does not yield the correct result. 
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Recently, Clarke and Pringle*^^ have shown that application of the mean free 
path theory to a rotating gas can indeed lead to the correct formula, Q- They 
consider the case in which the shear velocity of the gas is sufficiently smaller than 
the thermal velocity of the molecules. Their argument employs the inertial frame 
and approximates molecular orbits by straight lines. 

We have been inspired by Ref. 8) and, focussing on that paper, show here that 
their result can be derived with a simpler calculation procedure and a physically 
clearer picture. 



§2. Non-rotating, linear shear flow 



We consider, for simplicity, an inertial system in two dimensions, in which there 
exists a parallel flow in the y direction. We consider a point S in the flow and define 
its coordinates as (xo,0) (see Fig. 1 with U' < 0, where U' is the velocity gradient 
of the flow). An observer moves with S, and our argument employs his rest frame, 
which is also an inertial frame. As seen by the observer, the velocity of the flow can 
be written Au = (0, {x — xo)U'). Let us consider a point E that is at a distance A/2 
from S, where A is the mean free path of gas molecules.*-* Hereafter, we assume that 
the flight lengths between succesive collisions for all molecules are the same, i.e. A, 
for simplicity. 

The velocity of the flow at the point E as seen from S is directed in the y 
direction, and its magnitude can be written as 

Au = —-XU' cos a, (3) 

where a denotes the angle between the lines SE and SA. The flow velocity at A is 
also Au, as is clear from Fig. 1. The situation here is similar to that in Ref. 8). 

Gas molecules are assumed to be ejected from the point E isotropically when 
observed from the frame moving with E, and at constant speed c.^^ The flow velocity 
and molecular velocity added give v. Application of cosine theorem to the triangle 
formed by Au, c andv gives 



= i;^ + {AuY - 2vAucos (^-^ - a 

Following Ref. 8), we assume ^ {Au)'^, ignore higher-order terms, and obtain 

/ IXU' \ / IXU' \ 
V = c I 1 cos a sm a I = c I 1 sm 2a I . (4) 



2 c J \ 4 c 

Gas molecules ejected from the point E pass through the point S without a change 
in their velocities or directions of motion. Equation (4) shows that the velocity 
distribution of the gas molecules at S is no longer isotropic but exhibits an oval 
shape whose major axis is inclined 45° toward the x-axis. 



*■* Although Clarke and Pringle chose this distance as A, we consider S to be the midpoint of 
the trajectory of a molecule and take this distance to be A/2. 



Kinetic Theory Viscosity 3 

y 



i 

Au 


i 


A-cos(a)/2 




/ ^ 




A 


m / 


a y 


S 














/^a / 








\ E 











Fig. 1. Schematic diagram showing velocities of gas flow and a gas molecule. If ere, E is the emission 
point, S the observational point, c the thermal velocity of a molecule, Au the mean flow velocity 
at E and A, A the mean free path, and a the angle between SE and SA. In this figure V < 0. 



The x-y component of the viscous stress tensor axy is the net y momentum, 
carried by the gas molecules through a line of unit length along the y-axis per unit 
time. The x and y components of the velocity are Vx = v cos a and Vy = v sin a, 
respectively. The y momentum is mvy, and the mass flux is nvx, where m and n 
are mass of a molecule and the number density of molecules, respectively. Thus, 
axy = mn{vxVy) , where (• • •) represents an average taken over velocity space. 

Because we assume that the molecule velocity at E has a constant value c in all 
directions, the distribution function is non-zero only on the circle with radius c, and 
therefore it must have the form f(v) ~ 5{v — c), where 5 is the delta function and 
V is the radial coordinate in velocity space. Integration of — c) over v-a space 
yields 2ttc. Then, using v in Eq. @, our distribution function is found to be 



f{v,a) 



2'KC 



1 W 



■ sin 2a 



(5) 
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Substituting Vx and Vy into axy = mn{vxVy) , we obtain 

(7 xy = mn (^v'^ sin a COS a"^ = mn / / f{v,a)v'^ sin a cos advda. 

Jo Jo 

Fixing a and integrating the above over v from to oo gives 

mnc2 /-^^ / 1 AC/' \^ 

Cot = — : / 1 — 7 sin 2a I sin2Q:aa 

An In \ 4 c 







- — / 1 — 7 sin 2a sin2aaa = ——pcXU , 

An Jq \ 4 c J 16'^ ' 



where p = nm is the density of the gas. This equation, together with formula (1), 
gives the viscosity coefficient r] = 3pcX/16. The kinematic viscosity v (= rj/p) is, 
from the above formula, = 3cA/16. The usual heuristic procedure based on mean 
free path theory gives the coefficient 1/3 (for a three-dimensional treatment) or 1/2 
for a two-dimensional treatment, instead of 3/16. 

As is clear from the above calculation, if the velocity distribution at S is isotropic, 
the viscous stress there will become zero. Accordingly, the presence of viscosity 
requires anisotropy in the velocity distribution. In our calculation, the anisotropy 
in the velocity distribution of molecules at the point S has been derived assuming 
isotropy of the velocity distribution of molecules in the frame of the gas at the point 
E. This procedure resembles that in which, using the Chapman-Enskog expansion of 
the Boltzmann equation, one first assumes an isotropic Maxwellian distribution as 
the zeroth-order approximation in order to derive an anisotropic velocity distribution 
as the first-order approximation. 

We emphasize that this isotropy is only a calculational means for obtaining 
the anisotropic velocity distribution. If necessary, the velocity distribution obtained 
using the formula Q can be taken as that at the point E to calculate a second-order 
correction for that at the point S. 



§3. Circular flow 



3.1. Calculations in a rotating frame 

We next consider gas rotating about an axis with the velocity u = (0, i?l7(i?)), 
as seen from the inertial frame. For the moment, we ignore the effects of the Coriolis 
force and gravity and continue to assume that the molecule orbit is a straight line. 
In this case, the situation is again similar to that depicted in Fig. 1. The center of 
the rotation, O, is located far to the left (not shown in the figure). Let us denote 
the distance OE by Re and the angle SOE by (j). If A ^ R, then (p ^ X/R, and so 
cos (/) PS 1. We also have OE = OA^ i.e. Re ~ i? — ^A cos a. In the above derivation, 
we have ignored terms of second and higher orders of X/R. 

The y component of the fiow velocity at the point E as seen from the observer 
located at the point S differs from that given in the formula If the observer at 
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S moves with the inertial frame adopted in Ref. 8), it becomes 

AiU = REf2{RE) cos (t)-RQ^ Re^{Re) - RO ^ + RQ')\ cos a. (6) 

The X component of the flow velocity at E in the inertial frame is Re^^{Re) sm.(j) ~ 
0{\/R), which cannot be ignored in the present approximation. 

The y component of the velocity, however, when the observer at the point S 
moves with the rotating frame having an angular velocity of Q{R), becomes 

ArU = Re{0{Re) - n{R)) cos (t> ^ Re{^^{Re) - ^{R)) ~ -]^RQ'\ cos a. (7) 

This can also be expressed as A^u = AiU + l7Acosa/2.^) Indeed, for rigid body 
rotation, where Q is constant, this shows, as expected, that ArU becomes zero, 
while AiU is clearly non-zero. In the rotating frame, the x component of mean flow 
velocity can be ignored, because it is Re{^{Re) — f^{R))sin(l) = 0(A/i?)^. This 
is the advantage of employing the rotating frame rather than the inertial frame. 
However, the rotating frame has the disadvantage of the appearance of the Coriolis 
force. The Coriolis force, however, does not introduce a problem, as discussed below. 

Both the inertial and rotating frames can, with correct calculations, yield the 
same result. The important point is which provides a simpler calculation. The 
rotating frame obviously insures a simpler calculation, and for this reason we adopt 
this frame here. 

Comparison of the formula ((TJ with (jS)) clearly shows that there is simple re- 
placement of U' (for a plane shear flow in the inertial frame) by Rf2' (for a circular 
flow in the rotating frame) . It is then clear that the viscosity formula for the rotating 
flow has the same form as the formula 0. 

3.2. Calculations in the inertial frame 

Clarke and Pringle,®-* employing the inertial frame, after complicated calculations 
derived the result 

/ 1 \RQ' 

Vy = c sin all sin a cos a 

" V 2 c 

where their A is replaced by our A/2. Similarly, Vx can be derived easily: 

/ 1 \Rn' 

Vx = ccos all sin a cos a 

V 2 c 

No further calculations are presented here, because they are the same as those given 
above, except that U' has been replaced by Rf2' , as described above. The result is, 
naturally, the same as that obtained in the rotating frame. 

3.3. Effect of curvature of molecule orbits 

To this point, we have ignored the effect of the curvature of the molecule orbits. 
When in the rotating frame, one experiences the Coriolis force. The Coriolis force 
acts on a molecule, changing its trajectory from linear to circular. The radius of 
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this circle, which is the same as the Larmor radius of a charged particle moving in a 
magnetic field, is represented by c/l7. This radius is, for an accretion disk, believed 
to be approximately equal to the disk thickness, H. Under the condition A <C -ff, 
the trajectory of a gas molecule can therefore be regarded as linear. With the result 
based on the Boltzmann equation, we find that the Coriolis force does not change 
the form of the viscosity formula but, rather, produces anisotropy in the viscosity 
coefficients and suppresses the viscosity coefficient in the direction perpendicular to 
the rotating axis.^^~^) This effect does not appear in the present approximation, 
however. 

The effect of gravity can be studied similarly. We consider the rotating frame 
moving with the point S, which exhibits Keplerian motion. When considered on a 
sufficiently small scale about the point S, the rotating frame can be approximated 
by the Hill coordinates. It is known that a test particle moves in an elliptical (not 
circular) orbit whose major axis to minor axis ratio is 2:1.^) The curvature radius of 
the orbit is of the order of c/i7, which is the same as the above result. In summary, 
the molecule/particle orbit can also, within the range of our approximation, be 
approximated as linear, even when gravity is important. 

When A is large, the suppression of the viscosity coefficient in the rotational 
plane occurs for two reasons: a) shortening of the effective mean free path of gas 
molecules due to their curved orbits, and b) reduction of the asymmetry in the ve- 
locity distribution at S, because of the curvature of the particle orbit. This effect 
reduces {v-xVy) and, therefore, the effect of viscosity. In the extreme case, the veloc- 
ity distribution becomes symmetric about the ,T-axis as is seen from the numerical 
simulation of Narayan et al.^^ The above effects of a) and b) are of the order of 
{X/Hf. 

§4. Discussion 

We now return to the original question that led to our previous paper.^^ This 
question is set up as follows. In accretion disks, the angular momentum increases 
outwards, while the angular momenta of the molecules are conserved during their 
motion. Why, then, is angular momentum transported outward in opposition to the 
angular momentum gradient? 

Let us consider two adjacent annuli in a Keplerian rotating gas disk. The gas in 
the inner annulus has a larger angular velocity and smaller angular momentum than 
that in the outer annulus. Note that the velocity distributions of the molecules in 
both annuli are much larger than the velocity difference between the two annuli. Of 
the molecules in the inner annulus, many have much larger specific angular momenta 
than the average specific angular momentum of the outer annulus (while, of course, 
many more have smaller ones). 

We thus find that the above calculation shows that the velocity distribution at 
the point S comes to possess an oval shape whose major axis is inclined 45° toward 
the X-axis. This means that, of the molecules originally present in the inner annulus, 
those whose angular momenta are larger than the average angular momentum of the 
outer annulus (such molecules are ejected in the positive y direction) will preferen- 
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tially be transported from the inner to the outer annulus. The contribution from 
those molecules that have smaller angular momenta and have reached the outer an- 
nulus is small. Conversely, of the molecules originally present in the outer annulus, 
those having angular momenta smaller than the average angular momentum of the 
inner annulus will preferentially be transported to the inner annulus. Accordingly, 
angular momentum indeed does flow against its gradient. 
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